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Abstract. The HEGRA system of imaging atmospheric 
Cherenkov telescopes provides for specially selected 
classes of events an angular resolution of better than 3'. 
By comparing the measured angular distribution of TeV 
gamma rays from the Crab Nebula with the distribu- 
tion expected on the basis of Monte Carlo simulations, 
and with measurements of gamma rays from the point 
source Mrk 501, we conclude that the rms size of the VHE 
gamma-ray emission region in the Crab Nebula is less than 
1.5'. 

Key words: ISM: supernova remnants - ISM: individual 
objects: Crab Nebula - Gamma rays: observations 



1. Introduction 

Stereoscopic systems of imaging atmospheric Cherenkov 
telescopes (lACTs) such as the HEGRA lACT system 
(Daum et al. 1997 ) allow to reconstruct the directions 
and energies of TeV gamma-rays with high precision. The 
analysis techniques, the control of systematics, and the 
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understanding of the angular resolution function of the in- 
strument (Piihlhofer et al. 1997 , Aharonian et al. 1999b , 
Hofmann et al. 1999) has progressed to a level that one can 
start to study the characteristics of extended sources on 
a scale of a few arcminutes. These scales start to become 
interesting in disentangling the emission mechanisms for 
Galactic TeV gamma-ray sources such as the Crab Nebula 
or the pulsar PSR B1706-44. In this paper, we present, as 
a case study, an investigation of the size of the emission 
region of the Crab Nebula at TeV energies. 

The Crab Nebula is one of the best-studied objects in 
the sky, in all wavelength regimes. It has been established 
as a TeV gamma-ray source by the Whipple group, using 
the imaging atmospheric Cherenkov technique (Weekes 
et al. 1989| , Vacanti et al. 1991), and has been studied 
with many other Cherenkov telescopes. The precise spec- 
tral shape of gamma-ray emission from the Crab Neb- 
ula has been the subject of a number of recent publica- 
tions (Hillas et al. 1998| , Tanimori et al. 199S, Konopelko 
1999a). The spectrum is consistent with a power-law ex- 
tending from a few 100 GeV out to energies of 50 TeV 
and beyond. Contrary to observations in the X-ray and 
GeV gamma-ray regimes, the TeV gamma-ray emission 
does not show a pulsed component attributable to a di- 
rect contribution from the Crab Pulsar; pulsed emission 
is below 3% of the DC flux (Aharonian et al. 1999c, Bur- 
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dett et al. 1999). The commonly accepted model for VHE 
gamma-ray production in the Crab Nebula assumes elec- 
tron acceleration in the termination shock of the pulsar 
wind at a distance of about 0.1 pc (0.2') from the pulsar 
(see, e.g., Kennel & Coroniti (1984), De Jager & Harding 
( |1992D , Atoyan & Aharonian ( |1996| ), Aharonian & Atoyan 
(1998)). The electrons diffuse out into the Nebula and 
produce a characteristic two-component electromagnetic 
spectrum: synchrotron emission dominates at most ener- 
gies up to about 0.1 GeV, whereas the inverse Compton 
process generates higher-energy gamma-rays with energies 
from the GeV range up to 100 TeV and beyond. From the 
relative strength of the two components, values for the 
average mag netic field of 15 - 30 nT h ave be en derived 
(Hillas et al. |l994 De Jager & Harding |1992| , Atoyan & 
Aharonian |l996|) . 



The Crab Nebula represents an extended source of 
electromagnetic radiation. Since the electrons loose en- 
ergy as they expand out into the Nebula, primarily due 
to synchrotron losses, the effective source size is predicted 
to shrink with increasing energy of the radiation, with ra- 
dio emission extending up to and beyond the filaments 
visible in the optical, whereas hard X-rays and multi- 
TeV gamma-rays should be produced primarily in the di- 
rect vicinity of the shock (see, e.g.. Kennel & Coroniti 
( |1984|), De Jager & Hard ing ^1992| ), Atoyan & Aharo- 
nian ( |1996D , Amato et al. ( |1999 )). At TeV energies, a sec- 
ond production mechanism for gamma-rays could be the 
hadronic production by protons accelerated in the shock 



(Atoyan & Aharonian 1996) or resulting from decays of 



secondary neutrons produced in the pulsar magnetosphere 



(Bednarek & Protheroe 1997); gamma rays are produced 



in interactions with the surroun ding m aterial, e.g. in the 
filaments (Atoyan & Aharonian 1996 ). This contribution 
might be enhanced due to a trapping of protons in lo- 
cal magnetic fields associated with the filaments, increas- 
ing the interaction probability. Given that the size of the 
Crab Nebula - with its about 4' by 3' extension in the 
optical - is comparable to the angular resolution achieved 
for TeV gamma rays by the HEGRA system of imaging 
atmospheric Cherenkov telescopes (lACTs), a study of the 
size of the TeV emission region of the Crab Nebula is now 
possible with meaningful sensitivity. This paper reports 
such an analysis, based on data collected over the last 
years with the HEGRA lACT system. 

For comparison and later reference, we will briefiy sum- 
marize the existing information on the size of the Crab 
Nebula, as a function of the energy of the radiation. An 
obvious problem in such a compilation is that there is no 
unique definition of 'size'. For comparison with the TeV 
results given later, the most relevant quantity is an rms 
size, gained by approximating the intensity distribution 
by a two-dimensional Gaussian, or by directly calculating 
the rms width by projecting or slicing the intensity distri- 
bution along an axis, and averaging over direction s. Rm s 
width values based on 5 GHz radio data (Wilson |l972|) . 
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Fig. 1. Angular size of the Crab nebula at different fre- 
quencies. Full circles: rms size, averaged over directions. 
Open circles: half width at half maximum (HWHM, de- 
fined as FWHM/2), averaged over the long and short axis. 
Open squares: HWHM along the 300° direction. See text 
for references. The dashed line indicates the frequency de- 
pendence of the size of the (synchrotron-radiation) emis- 
sion region as given by Kennel & Coroniti (1984), and 
the full square shows the rms size predicted for inverse- 
Compton gamma-rays at TeV energy (Atoyan & Aharo- 
nian 1996). The dashed region indicates the rms size range 



of the filaments, the likely scale for hadronic production 
mechanisms. The triangles show the upper limits on the 
rms source size at TeV energies derived in this work. 



optical data from Woltjer (1957) as displayed in Wilson 
(|1972D , and 0.1 - 4.5 keV X-ray da ta (Ha rnden & Seward 
1984) are compiled in Hillas et al. ( 1998| ). Additional rms 



values were obtained for the 327 MHz radio data of Bi- 
etenholz et al. ( |1997| ) , and for the 22 - 64 keV X-ray data 
of Makishima et al. ( |l98lD . These data are shown in Fig. 
|l| as full circles. The bulk of the size values quoted in the 
literature refer to a different measure, the full width at 
half maximum (FWHM), which unfortunately in case of a 
structured intensity distribution depends also on the res- 
olution of the instrument. In some cases, data are only 
available along a specific direction, and do not allow to 
average over the long and short axis of the Nebula. In 
particular in earlier X-ray data, the contributions of the 
pulsar and of the surrounding nebula are not separated. 
Fig. |l| includes (as open circles) FWHM size data, aver- 
aged over the long and short axis, at radio wavelengths 
(Wilson |1972| ) and in the optical, NUV, FUV, and X-ray, 
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as given in Hennessy et al. (1992). Also shown are X-ray 
data compiled in Ku et al. (1976), which refer to the width 
along the 300° direction (open squares). A clear trend 
for decreasing source size with increasing energy is evi- 
dent, considering either the rms size values, the averaged 
FWHM size values, and the fixed-direction X-ray widths. 
Included as dashed line is the frequency dependence of 
the synchrotron emission region as sketched in Kennel & 



Coroniti (1984). The size of the emission region for inverse- 



Compton TeV gamma rays can be predicted using the 
average magnetic field to relate synchrotron photon ener- 
gies to electron energies and to inverse- Compton gamma 
rays; from such arguments, one concludes that the size 
for TeV gamma-rays should correspond to the X-ray size. 
The rms size predicted for the TeV gamma-ray emission 
region by the detailed calculations of Atoyan & Aharonian 
( |l996D is included in Fig. |[ Basically, inverse-Compton 
TeV gamma-rays should emerge from the toroidal X-ray 
emission region clearly visible in the ROSAT data (Hester 
et al. 1995| ) and in the recent Chandra image (Weisskopf 
et al., 2000), as already speculated earlier by Aschenbach 
& Brinkmann (1975). The projected semi-major and semi- 
minor axes of the emission torus are 38" and 18" , respec- 
tively. Due to the nonuniform strength of X-ray emission 
along the torus, the resulting emission profile is roughly 
elliptical, and its center is shifted N relative to the pulsar 
location by about 0.3'. Hadronic production mechanisms 
are expected to generate larger source sizes, of the scale 
of the size of the nebula (shaded region in Fig. |l|) . 

2. Observations of the Crab Nebula with the 
HEGRA CT system 

The HEGRA system of imaging atmospheric Cherenkov 
telescopes is located on the Canary Island of La Palma, 
on the site of the Observatorio del Roque de los Mucha- 
chos. The telecope system consists of five telescopes, with 
a mirror area of 8.5 m^ and a focal length of 5 m. A sixth 
prototype telescope is operated in stand-alone mode and 
is not used for the analysis presented here. The system 
telescopes are arranged at the corners and in the center 
of a square of about 100 m side length. The alt-azimuth 
mounted telescopes are equipped with cameras consist- 
ing of 271 photomultipliers (PMTs). Each PMT views an 
area of the sky of 0.25° diameter; the field of view of each 
camera is about 4.3°. Cherenkov images of air showers 
are recorded whenever two telescopes trigger simultane- 
ously; the trigger condition requires that two neighboring 
PMTs exhibit signals equivalent to 8 or more photoelec- 
trons. Typical trigger rates are around 15 Hz, for an en- 
ergy threshold of 500 GeV for vertical gamma rays. Details 
about the HEGRA lAC Ts an d their performanc e can b e 
found in Daum et al. ( 1997 ), Aharonian et al. ( 1999b| ), 



1996, initially with three telescopes, later with four and 
since late 1998 with the complete set of five telescopes. 
For this analysis, only data taken in the years 1997 and 
1998 were used, acquired with at least four telescopes un- 
der good weather conditions; in order to be able to com- 
pare with earlier Mrk 501 data taken with four telescopes, 
data from the fifth telescope was not used in the most re- 
cent five-telescope data sets. The quality-selected data set 
amounts to an integral observation time of about 155 h, 
and includes about 6.3 million events. For the final analy- 
sis, only data taken at zenith angles of less than 30° were 
included, with about 3.5 million events remaining. 

The Crab Nebula was observed in the so-called wobble 
mode, with the source offset by 0.5° in declination relative 
to the telescope axes. The sign of the offset alternated 
every 20 min. A region offset by the same amount, but in 
the opposite direction, is used for background estimates, 
avoiding the need for special off-source observations. Since 
the stereoscopic reconstruction of air showers provides an 
angular resolution of typically 6', signal and background 
regions are well separated. 

The techniques for data analysis are similar to those 



Hermann (|1995|), Buhan et al. ( |1998D . 

The Crab Nebula was observed in each season since 
the HEGRA lACT system commenced operation in late 



documented, e.g., in Aharonian et al. (19991 ,1999d). Di- 
rection and impact point of an air shower are recon- 
structed from the stereoscopic views of the shower. Based 
on the measured impact point and the known (Aharonian 
et al. 1999a) distribution of Cherenkov light as a function 
of the distance to the shower axis, an energy estimate is 
derived, with a typical resolution of 20%. Gamma-ray can- 
didates are selected on the basis of image shapes. Given 
the core distan ce and the intensity (the Hillas size pa- 
rameter (Hillas |l985 )), the expected width of a gamma 
ray image is determined. The measured width values are 
normalized to this value, and averaged over telescopes. A 
cut on the resulting mean scaled width < 1.2 retains most 
gamma-rays, but rejects the bulk of the cosmic-ray show- 
ers. 

In detail, the reconstruction of shower geometry differs 
somewhat from the techniques used so far. Whereas the 
normal reconstruction procedure combines images from all 
telescopes regardless of their quality, the new procedure 
assigns - on the b asis of the Monte Carlo simulations of 
Konopelko et al. (|l999b| ) - errors to the relevant image 
parameters (the location of the image centroid and the 
orientation of the image axis). These errors depend on the 
intensity and the shape of the images and are propagated 
through the geometrical reconstruction, resulting in error 
estimates (or, to be precise, a covariance matrix) for the 
shower parameters. Details of the algorithm are given in 
Hofmann et al. (199£). Depending on the characteristics 
of an event, an angular resolution between 2' and more 
then 10' is predicted, with average value slightly below 6' 
0. The ability to select subsets of events with better-than- 



^ Here and in the following, "angular resolution" is defined 
as the Gaussian width a of the distribution of reconstructed 
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average resolution will be used extensively in the analysis 
of the size of the VHE emission region in the Crab Nebula. 

3. The limit on the size of the emission region 

Even if the size at radio wavelengths - about 1.3' rms 
- is used as a most extreme possibility for the size at 
TeV energies, the source size as is still smaller than the 
angular resolution CTq for the best subsets - about 2' to 3'. 
Therefore, one cannot expect to generate a detailed map of 
the source. Instead, an extended emission region of (rms) 
size (Ts would primarily show up as a slight broadening of 
the angular distribution of gamma-rays, beyond the value 
determined by the experimental resolution: 



or, for (Ts small compared to ctq 



Aa 



2al 



(1) 



(2) 



For an intrinsic resolution of 3' (in a projection) and a 1.5' 
rms source size, one would find a 3.4' wide angular distri- 
bution; for the 6' resolution, the resulting width is 6.2'. In 
order to positively detect a finite source size, or to derive 
stringent upper limits, one has (a) to measure the width 
of the angular distribution of gamma-rays with sufficient 
statistical precision, and (b) to quantitatively understand 
the response function of the instrument at the same level, 
and to control systematic effects which influence the res- 
olution. 

For a given number n of events from the source, and 
ignoring for the moment the effect of background un- 
der the signal, the statistical error on the width of the 
(projected) angular distribution is Act/ct — l/(2n^/^). In 
case the angular distribution is consistent with a point 
source, Eqn. ^ then implies a (1 standard deviation) up- 
per limit Us < tJoTi"^/*. Therefore, one will want to min- 
imize (To by selecting a subset of events with particularly 
well-determined directions, even at the expense of event 
statistics. More critical are, in general, systematic errors. 
Pointing imperfections, changes in mirror alignment, etc. 
can cause differences in the angular resolution between 
data sets taken at different times, under different condi- 
tions, and between the data and the simulation. If, e.g., 
the intrinsic resolution of the instrument is known and re- 
producible to 10%, the minimum source size which can be 
reliably detected is 0.46cto, according to Eqs. 0or |[ Once 
more, a selection towards small (Tq is preferred. Among 
gamma-ray events, we find that about 1% of the events 
have a predicted resolution below 1.8' (0.03°), 6% below 
2.4' (0.04°), 15% below 3' (0.05°), and 60% below 6' (0.1°), 
respectively. In particular the samples with resolutions 
better than 2.4' or 3' combine good angular resolution 

shower directions, projected onto one axis of a local coordinate 
system. 
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Fig. 2. Angular distribution of reconstructed showers in 
a local coordinate system centered on the Crab pulsar. 
Events were selected on the basis of the predicted angular 
error (< 3' in both directions) and on the basis of image 
shapes {mean scaled width < 1.2). 



with acceptable statistics. The cuts on angular resolution 
have the additional benefit on enhancing the gamma-ray 
sample relative to the cosmic-ray background; cosmic-ray 
showers generate more diffuse images and have a worse 
angular resolution. Fig. 2 shows the angular distribution 
of events retained after a cut at 3' resolution in both pro- 
jections, applying only very loose additional cuts on event 
shapes (a cut on the mean scaled width at 1.2, which re- 
tains over 80% of the gamma-ray events). The selection 
also biases the sample towards higher energies, since high- 
energy events produce more intense images, with smaller 
errors on the image parameters. In the overall data sam- 
ple, the median energy of reconstructed events is 0.9 TeV; 
after a cut on the resolution at 3', this value rises to 2.0 
TeV. 

While the evaluation of statistical errors is straight 
forward, the control of systematic errors is more difficult. 
The pointing of the telescopes is referenced to and cor- 
rected offiine on the basis of star images (Piihlhofer et 
al. 1997), and the achievable pointing precision has been 
investigated in considerable detail. We are confident to 
achieve a pointing deviation of less than ±0.5' in each co- 
ordinate. Indeed, when the Crab data set was subdivided 
into seasonal subsets, the reconstructed source position 
was in all cases consistent with the nominal source loca- 
tion - the position of the Crab pulsar. 

To evaluate the level at which the angular resolution 
is understood, we use the sample of gamma-rays from 
the AGN Mrk 501 as a reference set (see Aharonian et 
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Table 1. Width of the angular distribution of events relative to the source, comparing the Mrk 501 and Crab data 
sets with Monte Carlo simulations using the measured gamma-ray energy spectrum as an input. Data sets are selected 
according to the estimate of the angular resolution, as provided by the shower reconstruction algorithm. The quoted 
width values are derived using a Gaussian fit to the projected angular distribution. For the last two rows of the table, 
only the central part of the distribution is fit; there are significant non-Gaussian tails (both in the Monte Carlo and 
in the data). 



Selection on 


Mrk 501 MC 


Mrk 501 data 


Crab MC 


Crab data 


angular resolution [arc min] 


[arc min] 


[arc min] 


[arc min] 


[arc min] 


< 2.4 


2.43 ±0.05 


2.46 ± 0.06 


2.41 ±0.05 


2.41 ±0.14 


< 3 


2.81 ±0.04 


2.83 ±0.05 


2.81 ±0.04 


2.70 ±0.10 


< 6 


3.58 ±0.03 


3.63 ± 0.04 


3.64 ± 0.04 


3.70 ± 0.09 


all events 


4.23 ± 0.04 


4.26 ± 0.04 


4.30 ± 0.04 


4.37 ±0.10 



al. ( 1999b , 1999d ) for details on this data set), assuming 
that Mrk 501 represents a point source. Table |l|, columns 
2,3 compare the measured angular distribution for differ- 
ent subsets of events with the Monte-Carlo predictions. 
'Angular resolution' again refers to the Gaussian width 
of the projected angular distribution of events. Excellent 
agreement between data and Monte-Carlo is seen for all 
data sets. We note that the resolution estimates given 
by the reconstruction algorithm are low by 10% to 15%, 
for the samples selected for good resolution. The < 2.4' 
sample, e.g., should show a 2.2' resolution, compared to 
the measured value of 2.46'. Given the relatively crude 
parametrization of image parameter errors used in the re- 
construction (Hofmann et al., 1999), a deviation at this 
level is not unexpected, and in any case the effect is fully 
reproduced by the simulations. 

After these preliminaries, we can now address the Crab 
data set. Fig. ||(a) shows the background-subtracted an- 
gular distribution of reconstructed showers relative to the 
source direction, projected onto the axes of a local coor- 
dinate system, after a cut on the estimated error of less 
than 3'. Superimposed is the corresponding distribution 
of Monte-Carlo events, generated with the measured Crab 
spectrum. Fig. |^(b) shows the corresponding comparison 
with gamma-rays from Mrk 501. In particular after the se- 
lection on good angular resolution, the two event samples 
are very similar in their characteristics (mean predicted 
resolution, mean number of telescopes contributing to the 
reconstruction, etc.) and can be compared directly, de- 
spite the differences in the energy spectrum of the two 
sources. Table 0, cols. 4,5 list the widths of the distri- 
butions for the Crab gamma-rays, and the corresponding 
simulations. In general, we find, within the statistical er- 
rors, good agreement between the Crab and Mrk 501 data 
sets, and between Crab data and Monte-Carlo. 

Another option to allow a direct comparison of the 
Mrk 501 and Crab data sets is not to look at the angu- 
lar deviations (x, y) between a gamma-ray and the source, 
but rather to normalize these quantities with the angular 
resolutions (Jo,x,'^o.y predicted event-by-event by the re- 



construction algorithm. Ideally, one expects to see a Gaus- 
sian distribution of unit width, independently of the en- 
ergy spectrum, the number of telescopes active for a given 
data set, etc. The observed distributions are indeed Gaus- 
sian; as mentioned above, their widths differ by up to 10% 
to 15% from unity, depending on the selection of events. 
Most importantly, however, these width values are consis- 
tent between all four sets of events (Mrk 501 MC, Mrk 
501 data, Crab MC, Crab data), within errors of 1.5% or 
less for the large-statistics sets where all events are in- 
cluded. This agreement demonstrates that there are no 
uncontrolled systematic effects between the two experi- 
mental data sets, or between the experimental data and 
the simulations. 

Since the width of the angular distribution of gamma- 
rays from the Crab Nebula is consistent with the expected 
width, and with the width observed for Mrk 501, we can 
only give an upper limit on the source size. Taking into ac- 
count the statistical errors on the Crab sample and on the 
reference samples, we find - following Caso et al. (1998) 
- 99% confidence level upper limits of 1.0' for the sample 
with a cut at 3' resolution, and 1.3' for the < 2.4' sam- 
ple. To be conservative, and since it was always planned 
to use the < 2.4' sample as a safest compromise between 
statistical and systematic uncertainties, we adopt the 1.3' 
limit. Adding in additional systematic uncertainties due 
to pointing precision, we quote a final limit of 1.5' for the 
rms source size at a median energy of 2 TeV. 

A possible contribution of gamma-rays from hadronic 
processes is generally expected to be most relevant at 
higher energies, in the 10 TeV to 100 TeV range. There- 
fore, the source size was also studied for an event sample 
with reconstructed energies above 5 TeV; data are con- 
sistent with a point source and the corresponding limit 
on the rms source size is 1.7'. Limited statistics prevent 
studies at even higher energies. 

The elliptical shape of the X-ray emission region sug- 
gests to perform the same analysis in a rotated coordinate 
system, with its axes aligned along the major and minor 
axes of the X-ray profile. Results obtained for the width 



6 



F. Aharonian et al.: Study of the VHE emission region in the Crab Nebula 



> 

o 
d 



ill 
> 

o 
6 



300 
250 
200 
150 
100 
50 


300 
250 
200 
150 
100 
50 






(a) 






T 


^^^^ V> -L CA KJ KA.CX.VCM. 


: i 


<> MC 




v 


r^TTff,Tiimf'T,TTi , , , , 1 , , , , 1 , , , , 1 , , 1 1 iSP^fTTsTT^ 



-25-20-15 -10 -5 5 10 15 20 25 

x,y [arc min] 



(b) 
Crab data 

Mrk 501 
data 



.t +. 



-25 -20-15 -10 -5 5 10 15 20 25 

x,y [arc min] 

Fig. 3. (a) Full points: background-subtracted angular 
distribution of gamma-rays, projected onto the axes {x, y) 
of a local coordinate system centered on the Crab pulsar. 
Events are selected to provide an angular resolution of 
better than 3'. Superimposed, as open points, the corre- 
sponding distribution of Monte-Carlo events, normalized 
to the same area, (b) Comparison of the angular distribu- 
tions of gamma rays from the Crab Nebula (full points) 
and from Mrk 501 (open points). 



in the major and minor direction do not show any signif- 
icant difference. Given the fact that the limits are large 
compared to the (rms-)size of the X-ray emission region, 
this observation is not surprising. The TeV source is re- 
constructed 0.2' from the location of the pulsar, and is 
consistent both with the location of the pulsar, and with 
the center of gravity of the X-ray emission region, within 
the systematic errors in the telescope pointing (less than 
±0.5' in each coordinate). While the statistical error alone 
is small enough to resolve a shift of 0.3' as observed in 



the X-ray image, current systematic errors prevent such a 
measurement at TeV energies. 

The limits obtained in this work are included in Fig. 

0- 

We note in passing that also the widths of the distribu- 
tions obtained for the AGN Mrk 501 are consistent with 
MC expectations, indicating the absence of a halo on the 
arcminute scale. Potential halo types for AGNs include 
wide pair halos (Aharonian, Coppi & Volk 1994) or nar- 
row halos caused by intergalactic magnetic fields (Plaga 
1995). For a relatively near source such as Mrk 501, a pair 
halo would be much wider than the field of view of the 
camera and could not be detected as an increased appar- 
ent source size; the other type of halo speculated in (Plaga 
1995) would be well below our resolution. We do not see 
any time-dependence of the source size, or any correla- 
tion with the TeV gamma-ray fiux. Details will be given 
elsewhere. 

4. Concluding remarks 

The size limits given in this work illustrate the precision 
which can nowadays be reached in TeV gamma-ray astro- 
physics; a number of potential galactic and extragalatic 
sources are predicted to be extended sources on this scale. 

A few remarks concerning the interpretation of the lim- 
its: the values quoted above refer to the rms source size, 
with the implicit assumption that source strength is dis- 
tributed over an area which is similar to, or small com- 
pared to the angular resolution of the instrument, such 
that the convolution of the source distribution and the 
Gaussian response function can again be approximated 
by a Gaussian distribution. This is obviously the case for 
a roughly Gaussian source, and was explicitly checked for 
two alternative distributions, namely (a) the case that the 
source strength is uniformly distributed over the surface 
of a sphere of radius r, resulting in an rms source width 
of r/\/3, and (b) the case that the source is uniformly dis- 
tributed inside a sphere of radius r, with an rms size of 
r/^/l. The limit on 1.5' rms size can be reliably translated 
into a limit on the radius of a surface source of 2.6', or on 
the radius of a volume source of 3.4'. We note, however, 
that one can construct source models with a large rms 
source width which would not be detected by our method. 
One example is a combination of a localized source at the 
center with a weak halo with an extension of a few degrees 
or more. Such a faint, diffuse halo could not be detected, 
yet formally results in a large rms width of the source. 
However, given that it is hard to imagine to have TeV 
emission even beyond the radio emission region, such sce- 
narios are hardly relevant for the Crab Nebula. 

The limit on the size of the TeV emission region of 
the Crab Nebula is, by a factor around 4, larger than the 
size predicted by the standard inverse Compton models 
for gamma-ray production in the nebula. The limits, how- 
ever, approach the sizes expected for hadronic production 
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models, where high-energy gamma-rays are produced by 
nucleon interactions, more or less uniformly throughout 
the nebula. 
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